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ABSTRACT 

We have studied the molecular hydrogen energetics of the edge-on spiral galaxy NGC 891, using 
a 34-position map in the lowest three pure rotational II2 lines observed with the Spitzer Infrared 
Spectrograph. The S(0), S(l), and S(2) lines are bright with an extinction corrected total luminosity 
of ^ 2.8 X 10^ Lq, or 0.09% of the total-infrared luminosity of NGC 891. The II2 line ratios are nearly 
constant along the plane of the galaxy - we do not observe the previously reported strong drop-off 
in the S(1)/S(0) line intensity ratio in the outer regions of the galaxy, so we find no evidence for the 
very massive cold CO-free molecular clouds invoked to explain the past observations. The H2 level 
excitation temperatures increase monotonically indicating more than one component to the emitting 
gas. More than 99% of the mass is in the lowest excitation (Te^; ~125 K) "warm" component. In the 
inner galaxy, the warm II2 emitting gas is ^15% of the CO(l-0)-traced cool molecular gas, while in the 
outer regions the fraction is twice as high. This large mass of warm gas is heated by a combination of 
the far-UV photons from stars in photo-dissociation regions (PDRs) and the dissipation of turbulent 
kinetic energy. Including the observed far-infrared [01] and [CII] fine-structure line emission and 
far-infrared continuum emission in a self-consistent manner to constrain the PDR models, we find 
essentially all of the S(0) and most (70%) of the S(l) line arises from low excitation PDRs, while most 
(80%) of the S(2) and the remainder of the S(l) line emission arises from low velocity microturbulent 
dissipation. 

Subject headings: galaxies: individual (NGC 891) — dust, extinction — infrared: galaxies — infrared: 
ISM — ISM: molecules 



1. INTRODUCTION 

Half the interstellar medium in late type spiral galax- 
ies is in molecular gas clouds, and about 90% by number 
of the atoms within these clouds are tied up in the H2 
molecule. Stars form within molecular clouds, and the 
pure rotational lines of H2 are important coolants en- 
abling cloud collapse. The cooling radiation from these 
lines is also important to the formation of the giant 
molecular clouds (GMCs) themselves. Despite its abun- 
dance and importance, very little is known about cool 
H2 on galactic scales. This is because H2 has no dipole 
moment, so that dipole radiation from its low-lying ro- 
tational energy levels is forbidden. The low lying lev- 
els instead radiatively decay by emitting relatively weak 
quadrupole (AJ — 2) radiation. Therefore the pure ro- 
tational lines of H2 are challenging to observe. Three of 
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the low lying lines: the S(l) (IT.O^m), S(2) (12.3^m), 
and the S(4) (8.0/im) lines are transmitted through the 
Earth's atmosphere, and have been detected from r egions 
as diverse as supernovae remnants (jRichter et al.l[T995.) . 
photo-dissociation regio ns (PDRs) associate d with OB 
star formation sites (e.g. iParmar et al.lll99lD. and disks 
enveloping young stellar objects (e.g. lBitner et al.| [2007'). 
However, detection of the lowest lying S(0) (28.2/im) line, 
and detection of these lines from cool molecular clouds in 
external galaxies awaited the advent of the Short Wave- 
length Spectrometer s (SWS) on the Infrared Space Ob- 
servatory (ISO, e.g. [We^Hn^^nder^rKllll). A 
few years later, the Infrared Spectrograph (IRyj) of the 
Spitzer Space Telescope became available and the detec- 
tion of these lines from a wide variety of external galax- 
ies became routine fe.g. lArmus et all 120041 : iSmith et al.l 
I2OOI IBernard-Salas et al.ll2009n . 

H2 exists in two varieties: para (nuclear spins anti- 
aligned) and ortho (spins aligned). The relative ortho- 
to-para {o/p) abundances refiect the local gas tem- 
perature since proton exchange reactions with H"*", H 
and H j" are effective in thermalizing the o/p ratio 
(jSternberg fc Neufeldlll999l and references therein). At 
high temperatures (T>200 K) the o/p ratio is set by 
the ratio of statistical weights at 3:1. Since there is no 
radiative coupling between the ortho and para species, 
when analyzing radiative transitions they can be treated 

^ The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. Sup- 
port for this work was provided by NASA through Contract Num- 
ber 1257184 issued by JPL/Caltech. 



as independent molecules. The lowest lying transitions 
for the para species are: J=2-0 (S(0), 28.2/xm), and 4-2 
(S(2), 12.3/xm), and for the ortho species: J=3-l (S(l), 
17.0/im), and J=5-3 (S(3), 9.7/iin). With their very small 
Einstein A coefficients the H2 level populations are easily 
thermalized at molecular cloud densities, and the emitted 
lines are nearly always optically thin. Therefore, the line 
intensities trace the column density of the emitting level, 
and since the line emitting levels are thermalized and 
lie hundreds of degrees above the ground state (larger 
than typical molecular gas temperatures), the line ra- 
tios within a species should reflect the gas temperature, 
Tgas- For example, the S(2)/S(0) line ratio gets larger 
by a factor of 150 as Tgas goes from 70 to 100 K. The 
para J=2 level lies just 510 K above the ground so it can 
be emitted by cool {Tgas > 70 K) molecular gas and is 
an important coolant for molecular clouds in, for exam- 
ple, the early stages of collapse into protostars (T^lOO 
K, UHa -lOOcm-^). 

Are the rotational lines of H2 detectable from nor- 
mal star forming galaxies, and if so, is the emit- 
ting gas important? Work with the ISO-SWS demon- 
strated that the answer to both questions is an em- 
phatic yes. Large scale S(0) and S(l) line emis- 
sion is reported both from the face-on Sc galaxy 
NGC 6946 (Vale ntiin et all [TOOl. and the edge o n 
Sb galaxy NGC 891 (IValentiju & van derWerJ IllH. 
The strong emission indicates surprisingly large masses 
of warm molecular gas. Spitzer IRS studies con- 
firm that strong H2 rotational line emission and large 
masses of warm molecular gas are common in a 
wide var iety of galax ies including nucl ei of starburst 
galaxies (iDevost et all ^ 2004: Bcrnard-Sa las et al.l [2009( ) 
and ULIRGs ( iHigdon ct al. 2006b), in disk of galax- 
ies (l Roussel"et al. 20071. extragalactic regions such as 
Stephan's Quin tet (lAppleton et al.ll2006[ ). as well as tidal 
dwarf galaxies (jHigdon et al.ll200 6a'). 

NGC 891 is a n earby (d = 9.5 Mpc 
Ivan der Kruit fc Searle I (|1981| )) Sb galaxy that is 
presented to us nearly edge-on and is thought to be 
a close Milky Way analogue. Since an edge-on galaxy 
presents the largest column densities to our telescopes, 
and the H2 rotational lines are optically thin, NGC 891 
provides the best opportunity to study warm H2 in 
the outer reaches of ga l axies. With the ISO-SWS, 
I Valentiin fc van der Werj ([1999) detected the S(0) and 
S(l) lines at eight positions along the plane of NGC 891. 
The lines were surprisingly strong even in the outer 
regions of the stellar disk. Assuming o/p = 3 they find 
warm {Tgas '~130 to 230 K) molecular cloud component 
pervades the galactic disk with a total mass comparable 
to both the atomic HI and the cold molecular gas as 
traced through its CO emission: M^arm h^ ^^coid co 
~ Mhi- It would be a surprise if half the molecular ISM 
in the galaxy is warm, and it is a challenge to provide 
a heating mechanism for this gas. Furthermore, based 
on a drop in the S(1) /S(0) line intensity ratio a t radi i 
greater than 8 kpc, I Valentiin fc van der Werj ()1999f ) 
argue for an additional cool (Tgas < 90 K) component 
which becomes dominant at large radii. With o/p ~2 to 
3, the cooler component is massive: Mcooi H2 '^ 5-15 x 
Mhi in the outer disk. The cool molecular gas mass is 
so large that it can be a substantial fraction of the " dark 



matter" required to drive the rotation curve out to the 
edge of the opitical disk of NGC 891. This is particularly 
intriguing due to theoretical work that suggests a large 
fraction of the "dark" baryonic m atter in galaxies may 
be in CO free H2 clouds (e.g. iPfenniger &: CombesI 
|1994[ ). The detection of a massive cooler component, if 
proven correct, is a very significant result. These results 
need verification so we have repeated the observations 
of the S(0) and S(l) lines along the plane of NGC 891 
using the more sensitive spectrograph on the Spitzer 
Space Telescope (Hpuck et al. 2004). In addition, the 
important S(2) line is included in the large instantaneous 
bandwidth of the Spitzer IRS. The addition of the S(2) 
line mapping provides a key element in tracing the gas 
excitation. 

2. OBSERVATIONS 

The galaxy was observed on 2004 August 7 as part of 
the IRSguaranteed time (GT) program (PID 97) with the 
two high resolution (R~600) modules, Short-High (SH) 
and Long- High (LH). We also obtained supplemental, 
deep integrations in selected positions through a Spitzer 
Cycle 4 open time (OT) program with PID 40877. For 
all observations, the SH slits covers the 9.7-19.5/im range 
and it is 4.7" x 11.3" in size, while LH covers the 18.8- 
37.2yLtm range and its size is 11.1" x 22.3". The width of 
each slit was designed so that it is equal to the first Airy 
ring of an unresolved point s ource at their short wave- 
length fsee lHouck et al.l l2004'). We obtained spectra for 
a total of 27 positions symmetrically sampled around the 
nucleus of the galaxy at a uniform spacing of 30" (1.53 
kpc) with a position angle of 23 degree along the galactic 
plane of NGC 891. Six positions off the galactic plane to 
the northwest were also included to search for possible 
extended H2 emission. For each pointing, the on-source 
integration time was 60 and 28 sees for SH and LH re- 
spectively. In addition two off positions to the northeast 
and northwest away from the galaxy, which were used for 
background subtraction, were also obtained with twice 
the integration time. The IRS pointings are indicated in 
Figure 1 with circles of 25" in diameter, slightly larger 
than the length of the LH slit. The enumeration of the 
pointings follows the notation of Table 1, which provides 
the RA and Dec of each pointing along with its distance 
in arcsec along the disk of the galaxy from the central 
pointing (ID=:1). Within the OT program, we obtained 
on-source integration times of 750 and 900 seconds for 9 
positions (positions 1,6,8,10,12,20,22,24 and 26 in the no- 
tation of Figure 1) and 5 positions (positions 1,6,8,20,22) 
in SH and LH respectively. We also re-observed the off 
possitions for background subtraction. 

Given the fact that IRS has no moving parts and that 
the long axis of the SH slit is perpendicular to that of the 
LH slit, we could not modify at will the orientation of the 
IRS slits on the sky. The slit orientation is determined 
by the ecliptic coordinates of the science target and the 
time of observations. By properly selecting the time of 
the observations, we were able to ensure that for the SH 
observations, which contain the S(l) and S(2) lines, the 
SH slit was placed with the long axis perpendicular to the 
plane of the galaxy. The LH observations, which trace 
the S(0) line, were obtained at the same time, with the 
LH slit placed with the long axis along the plane of the 
galaxy, since the two slits are perpendicular. 



The data were processed by the Spitzer Science Center 
(SSC) pipehne version 17.2. As it is discussed in detail in 
the IRS manual, for each pointing, two observations are 
made by placing the target on two nod positions along 
the slit. The 128 x 128 pixel detectors are read with A^ — 1 
successive non destructive reads before the final n-th read 
which also resets the array. For our observations N=16 
for SH and 8 for LH. The SSC pipeline starts with the 
128x128 X N SH data cubes and converts them to a 
128x128 two-dimensional image after linearization cor- 
rection, subtraction of darks, cosmic-ray removal, stray 
light and fiat field correction. Each image, called basic 
calibrated data (BCD) contains the ten echelle orders of 
the SH and LH slits respectively. Our analysis started 
from these BCD files. At each pointing, we subtracted 
the dedicated off-source background image in order to 
remove the sky emission from our slit. Then, the two 
dimensional images were cleaned with the IRSCLEAN 
package to remove the bad pixels as well as to correct 
the so called "rogue" pixels. These are pixels whose val- 
ues depend not only on the photons they collect, but 
also on the total number of photons hitting the whole 
128x128 array, in a nonlinear manner Ej Finally, the 
cleaned images were used to extract the spectra using 
the full slit extra ction method of th e IRS data analysis 
package SMART (jHigdon et al.ll2004[ ). The outputs from 
SMART produced one spectrum at each nod position, 
and thene the two nod position spectra are combined. 
The data from the edge pixels within each echelle order 
were manually cleaned by truncating a few pixels at the 
beginning and end of each order and when there was in- 
consistency, the blue end of the order was preferred. No 
scaling was needed between the adjacent grating orders 
within the same module. The final errors include both 
statistical and systematic effects, and we report them in 
Table 1. The overall absolute spectrophotometric uncer- 
tainty of IRS is 5%, so that spectrophotometric uncer- 
tainty is small compared with the signal-to-noise ratio of 
our spectra. 

3. RESULTS 

3.1. Line Morphology 

Table 1 lists the observed line fluxes and upper lim- 
its at all positions, while the locations on the galaxy 
where these spectra were obtained superposed on an 8/^m 
Spitzer/IRAC image of the galaxy, are indicated in Fig- 
ure 1. We also plot in Figure 2 the spectra around the 
S(0), S(l), and (S3) H2 lines from locations where at least 
one line has been detected. Both the S(0) and S(l) lines 
are clearly detected at 17 of the positions along the galac- 
tic plane. Figure 3a plots the line intensity distributions 
along the plane of galaxy. The lines have fairly uniform 
intensities in the inner galaxy, with some enhancement 
at the nucleus, beginning to fall off at 6.1 kpc NE and 
4.6 kpc SW of the nucleus to values below our detection 
threshold at regions beyond 13.8 kpc NE and SW of the 
nucleus. The S(2) line was only detected in 11 positions, 
all within 10 kpc of the nucleus, but this is very likely 
a sensitivity issue. The detection limit is similar for all 
three lines, but the S(2) is typically 4 times fainter than 
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the S(0) line and 1.7 times fainter than the S(l) line. 
Therefore, its distribution could well be nearly identical 
to the lower J lines. The line intensity ratios are also 
fairly constant across the galaxy, with the exception of 
enhanced S(l) line emission at the nucleus. 

We do not detect H2 line emission off the plane of the 
galaxy. However, our integration time is r ather mod- 
est, so perhaps this is not too surprising: iRand et al.l 
(2008) obtained substantially deeper integrations (nearly 
6.5 hours, or 400 times longer than ours) at two positions 
1 kpc off the disk, where th ey clearly d e tected the S(l) 
line. Based on these results IRand et al.l (J20085 estimate 
the scale height of the S(l) line emission is ~8.4 to 9.9" 
(390-490 pc, based on a distance of 9.5 Mpc). Correcting 
their results for the expected extinction of the lines by 
dust, they estimate that the intrinsic width of the warm 
molecular gas is -4.8 to 5.7" (220 to 260 pc). 

3.2. Comparison with ISO Results: No Evidence for a 
Massive Cold Component 

It is important to compare our Spitzer IRS re- 
sults to th ose from the ISO SWS study reported in 
(iValentijn et al. 1996). Howevever both the ISO and 
Spitzer maps are sparsely sampled, and neither was ob- 
tained with a large enough beam to detect all of the 
ex tra-planar emission if it is as extensive as the results 
of IRand et al.l ()2008D suggest. We first correct for the 
missing extra-planar emission, by assuming the appar- 
ent (exponent ial) scale h eight for the H2 line emission is 
the average of lRand et al. (2008) 's values: 9.1". We took 
Spitzer LH (S(0) line) data with the slit width extending 
±5.55" perpendicular to the plane, and SH (SI) and S(2) 
line data with the long axis of the slit expending ±5.65" 
perpendicular to the plane of the galaxy. For these cases, 
we calculate that 45% and 46% of the line fiux (inte- 
grated perpendicular to the plane) is detected within our 
LH and SH beams respectively. The (rectangular) ISO 
beams were 20" x 27" and 14" x 20" for the S(0) and S(l) 
lines respectively, and these beams were tilted so that 
the l ong axis of the rectan gle traced the plane of NGC 
891 ( Valentijn et al.l [19961 ). For these cases, and again 
comparing against an assumed scale height of 9.1" and 
integrated perpendicular to the plane, we calculate that 
45% and 53% of the line fiux was detected by the ISO 
SWS spectrometer. We apply these extra-planar correc- 
tions to both data sets for proper comparisons. We also 
must correct for flux not detected along the plane of the 
galaxy, since neither ours, nor the ISO data sets are fully 
sampled. The ISO SWS slit along the plane of the galaxy 
(point-to-point, along the diagonal) was 34" and 24" for 
the S(0) and S(l) lines respectively. Since the Spitzer 
data are only sampled every 30", we must interpolate 
between Spizer samplings to properly compare with the 
ISO data. For both the S(0) and S(l) data sets, we in- 
terpolate linearly between the Spitzer sampled regions, 
then integrate within a region covering the angular ex- 
tent of the ISO beam along the plane. This should be 
a good estimator for most of the galaxy, especially for 
the S(0) line, where the ISO and Spitzer beam sizes are 
nearly the same. A possible exception is the nucleus in 
the S(l) line, where the relatively much smaller Spitzer 
beam may emphasize any strong central peak. Figure 3b 
illustrates the comparison between Spitzer and ISO ob- 
servations including all the effects described above. For 



the inner (r < 7 kpc) galaxy our Spitzer line fluxes are 
about a factor 1.55 and 1.85 larger than those of ISO 
in the S(0) and S(l) lines respectively, but the Spitzer 
S(1)/S(0) line ratio is not much different than the ISO 
S(1)/S(0) ratio. Therefore, we derive similar excitation 
temperatures for the inner galaxy as those presented in 
iValentiin et al.l (|1996[ )fsee section 4.1.1 below). 

However, in the outer galaxy (r > 7 kpc) the S(1)/S(0) 
line ratios diverge. The Spitzer S(l) line flirx does not 
fall off nearly as rapidly as that measured by ISO, so 
that the Spitzer/ISO ratio grows to ^ 2.57 in the outer 
regions. In addition, the Spitzer S(0) line flux falls off 
more rapidly in the outer galaxy than does the ISO S(0) 
line flux, so that Spitzer/ISO ratio gets smaller, to ^ 
0.81 in the outer regions. Therefore, unlike our Spitzer 
results, where the S(1)/S(0) line ratio is nearly constant 
over the disk, the ISO results show a marked decrease 
in the S(1)/S(0) ratio for the 3 regions sampled at r > 
7 kpc from the nucleus. It is this change in line ratios 
for th e outer regions that led IValentiin fc van der Werj 
(|1999D to invoke a second, cool gas contributor to the 
observed S(0) line emission. Beam size is not the source 
of this discrepancy: since the ISO beam has its diago- 
nal along the plane of the galaxy, it was equally sensitive 
to extraplaner emission as the Spitzer beam. Further- 
more, it is not just the added S(0) line flux in the outer 
disk, but it is also the dearth of S(l) line radiation in the 
outer disk reported in the ISO observations that changes 
the observed S(1)/S(0) line ratio, hence the derived gas 
excitation. The nearly constant line ratios revealed by 
the Spitzer observations are consistent with warm molec- 
ular gas, with near uniform excitation conditions over 
the galaxy - there is no evidence for an additional cool 
gas component. Furthermore, as we will show below, 
this warm molecular gas component is far less massive 
(<13%) than the cold II2 emit ting molecular component 
> 2 2 X 10^ Mq) invoked bv IValentiin fc van der Werj 

nn). 

3.3. Extinction Correction 

To begin an analysis of mid-IR line emission from the 
edge-on galaxy NGC891, we need to first make correc- 
tions for the non-negligible extinction in the lines due 
to dust along the line of sight. The H2 line emission is 
likely to be highly correlated with that of the CO(l-O) 
line which is tightly confined to the plane of the galaxy, so 
that to make a good extinction correction, we need high 
spatial resolution proxies. Fortunately, there exist high 
resolution images of NGC 891 in the two primary tracers 
of the neutral ISM - the HI 21cm line which traces atomic 
hydrogen, and the CO(l-O) line which traces cold molec- 
ular hydrogen. We calculate the total column density in 
the molecular gas and atomic gas components within our 
Spitzer beams as a function of distance along the plane 
of the galaxy as follows. 

For the CO(l-O) d istribution, we use the data con- 
tained in Figure 2 of iScoville et al.l ()1993f ). These data 
were obtained with a 2.3" beam, but are plotted in their 
Figure 2 by integrating the emission perpendicular to the 
plane of the galaxy, then refering it to a 1" strip along the 
plane. To compare with our data, we therefore re-expand 
their strip to follow an exponential fall-off perpendicular 
to the plane of the galaxy with t he local CO(l - O) sca le 
height (also given in Figure 2 of IScoville et al.l (|1993l) ). 



The observed scale height (given as fuU-width-at-half- 
maximum) varies and corresponds to a 1/e scale height 
ranging from about 2.9" at the nucleus to ~ 4.3" in the 
outer regions. For regions just exterior (further removed 
from the nucl eus along t he g alactic plane) to those ex- 
plored bv iScovilleet al.l (|T993) (r>200") we have added 
in the IRAM 30 m te l escope (20" beam) observations 
of iGarcia-Burillo et al.l (|1992D . scaled assuming the line 
emission follows the 4.3" scal e heig ht observed in the 
outer galaxy bv iScovilleet al.l () 19931 ). We used a C0(1- 
0) integrated intensity to H2 column density conversion 
factor of X = 2 x 10^° cm~^/(K k m s~^) as for the solar 
neighborhood (jDame et al.ll200W^ . and derived an ex- 
pected hydrogenic column density within 1" strips par- 
allel to the galactic plane with an exponential fall-off as 
given in Scovil le et al . (1993). 

The HI column den sity is taken from the 21 cm line 
map of IRupenI (|1991f ) with a 7" beam. Since the scale 
height of the HI em ission is relatively large (~ 6", 
lOosterloo et al.l (|2007[ ). we estimate the HI column den- 
sity within our Spizer beams by direct inspection of the 
HI maps. Note that since the HI column is a small frac- 
tion of the H2 column (as traced by CO) in the the re- 
gions of the galaxy where extinction is important (within 
^ 5 kpc of the nucleus, see Figure 4a), a more detailed 
model of the HI distribution is not warranted. 

The total neutral gas mass is nearly equally divided 
into molecular (- 3.8 x lO'^ M©, Scoville et al. 1993), 
using X = 2 X 10^" cm~^/(K km s~^) , and atomic (^ 
4.1 X 10^ M0 lOosterloo et all I2007D components. We 
plot the calculated column density of hydrogen nuclei 
derived from the CO and HI observations, smoothed to 
our S(0) beam size along the plane of the galaxy in Figure 
4a. The column plotted only includes mass within our 
Spitzer beam which extends ± 5.55" perpendicular to the 
galactic plane: the mass missed at high galactic latitude 
is not included in this plot. The total column density 
of hydrogen nuclei within our beam peaks at 1.28 x lO^'^ 



cm at the nucleus and is above 2 x 10 



from -10 



to +10 kpc along the galactic plane. 

Such large columns lead to sizable extinction even for 
the S(0) line at 28/im. To find the extinction, we take the 
column density of hydrogen nuclie as traced by the CO 
distribution (modeled as an exponential fall-off) and the 
observed HI distribution within our beam to calculate 
the expected extinction as a function of distance from 
the galactic plane. The extinction calculation is broken 
up into 1" strips parallel to the plane. For each strip, we 
first find Aa, and then we calculate the local extinction 
correction, C(A) as defined below. We then average these 
extinction corrections within our beam to get an estimate 
for the total correction that must be apphed to the ob- 
served H2 line emission. Exti nction i s estim ated using 
the "Galactic" dust model of iDraind (|2003D for which 
A(28Aim)/N(H) ~ 8 x lO'^^ magcm^, A(17/im)/N(H) - 
A(12^m)/N(H) ^ 1.43 X 10"^^ magcm^, and make the 



^^ A different choice of X will not cfiange our results in any 
important way. For example, if X = 3 X 10 cm~^/(K km s~^), 
then the extinction corrected luminosities of the H2 lines go up 
10%, and the overall calculated warm (H2 line emitting) to cold 
(CO(l-O) line emitting) molecular mass fraction goes down to 11% 
from 16% (see section 4.2.2). The derived physical parameters of 
the gas, and its excitation mechanisms do not change significantly. 



correction for the hi gher gas to dust mass ratio observed 
in NGC 891 (260, see lAlton e t al. 200( | as compared with 
the Galaxy (160, see lSodroski et al.iil994D . We assumed 
a "mixed" extinction scenario, where the dust is mixed 
with the warm H2 gas, and for which the extinction cor- 



rection, C(A) is given by C(A) = t\{1 - e 



At the 



nucleus we derive a (beam averaged) extinction correc- 
tion for the S(0) line of C(28^m)-1.38, and for the S(l) 
and S(2) lines C(12/Ltm), C(17;Um)-1.71. The extinction 
correction as a function of offset from the nucleus is plot- 
ted in Figure 4b. The correction is significant for much of 
the inner (r j 5 Kpc) galaxy, and larger for the S(l) and 
S(2) lines, so that the extinction correction will change 
the line ratios. 

3.4. Line Luminosity 

The line fluxes change only gradually over the plane of 
the galaxy, and the slits likely enclose nearly all of the 
flux perpendicular to the plane. Therefore, to find the 
total line luminosity, we simply integrate the line fluxes 
along the plane of the galaxy, interpolating in a linear 
manner between samplings. For both the S(0) and S(l) 
lines, this amounts to integrating the observations and 
multiplying by the ratio of the sampling interval to the 
beam size along the plane: 30/22.3 = 1.35 and 30/4.7 = 
6.38 for the S(0) and S(l) lines respectively. The total 
line fluxes so obtained are F(S(0)) - 1.95 x 10"^^ W m"^ 
and F(S(1)) - 3.20 x lO^^^ W m^^. At a distance of 9.5 
Mpc, the corresponding line luminosities are 5.44 x 10^, 
and 8.94 x lO'^ Lq for the S(0) and S(l) lines respec- 
tively. The S(2) line is not as well sampled as the S(l) 
line due to sensitivity limits, but assuming the S(2)/S(l) 
line ratio for the detected positions (~ 0.59) is the same 
throughout the galaxy, the S(2) line luminosity would 
be ~ 5.27 x 10^ Lq. The total observed (and interpo- 
lated) power for these three low-lying lines is therefore 
- 1.97 X 10^ Lq, or 6.3 x 10"^ of the total IR (TIR, 3- 
1100 /im) luminosity, 3. 14 x 1 0^° Lp^ (flux densities from 
Brauhcr, Dale & and H eloul (pOOsf l. TIR as deflned in 
Dale fc Helou 1 (|2002i )). The observed low J molecular 
hydrogen lines are therefore ^ 1.5 times as bright rela- 
tive to the TIR continuum when compar ed to the s ample 
of 57 normal galaxies surveyed by Roussel et al.l ()2007f ) 
who find a typical ratio near 4 x 10~^ with a scatter of 
about a factor of two. 

Correction for extinction before integrating along the 
plane (factors of 1.16, 1.32, and 1.32 for the S(0), 
S(l) and S(2) lines respectively, section 3.3, above), 
which is bound to be more important for the edge-on 
NG C 891 than for the t ypically less inclined galaxies in 
the lRoussel et al] ()2007[ ) sample, the total H2 luminosity 
is ~ 2.51 x 10^ Lq or 8.0 x 10"'' of the TIR luminosity, o r 
~ twice the typical values found in lRoussel et al.l (120071). 

The question of line flux at high scale heights missed by 
the Spitzer beam is more difScult to quantify. If the scale 
height derived bv lRand et al.l (J2008D applies the flux we 
observe is roughly equal to that we miss at scale heights 
greater than that of our beam {^ 5.6"). The missed (high 
z) emission summed over the three lines is in this case 
therefore ^ 1.97 x 10^ Lq. The high-z emission would 
suffer very little extinction, so the the total luminosity 
in ah three lines would be ^ 2.51 4- 1.97 = 4.48 x lO'^ 
Lq in the case of a 9.1" (|Rand et al.ll2008D scale height. 



Another estimate is presented by presuming the H2 scale 
height is the same as that of the CO(l-O) line. Since 
the CO scale height is observed along the entire galactic 
plane, this is the correction we prefer. For this case, the 
flux exterior to our beam (without extinction correction) 
is just 16% of the flux within our beam (before extinction 
correction) so that the total luminosity in all three lines is 
~ 1.97x0.16-^2.51 = 2.83 X 10^ Lq. Correcting for both 
extinction and extra-planar emission, the H2 luminosity 
is therefore - 0.9 to 1.43 x 10"^ of the TIR luminosity. 

4. DISCUSSION 

4.1. Intensity Ratios 

To derive the physical properties of the emitting gas, 
we need to compute ratios between lines measured within 
different apertures. Inspecting Figure 3, it is clear that 
for the most part, the flux in the three lines changes 
gradually at the scale of LH slit length (22.3") along the 
plane of the galaxy. The S(l) and S(2) lines are observed 
with the long (11.3") axis of the SH slit perpendicular to 
the plane so that it matches the z-sampling of the S(0) 
beam. Therefore, to take ratios among the H2 lines, we 
need only correct for the beam size along the galactic 
plane. We do this by linearly interpolating between S(l) 
and S(2) samples (taken with a 4.7" beam along the plane 
of the galaxy), and then estimating the average intensity 
expected in the 22.3" S(0) beam from this interpolation. 
To obtain the physical parameters of the gas, we use the 
extinction corrected values of the fluxes contained within 
our Spitzer beams presented in Figure 5a. The extinction 
corrected line ratios are presented in Figure 5b. 

4.2. Gas Excitation: Analytical Models 
4.2.1. The Ortho to Para ratio 

The low lying rotational lines of H2 are optically thin 
and easily thermalized so that the lines and their ratios 
probe the temperature and column density of the emit- 
ting gas. The S(2) and S(0) lines both arise from the 
para species of H2 so that it in the high dens i ty lirn it 
(uH^ >nc„t(J=4) ~ 10^ cm- MLe Bourlot et al.l (|1999D ). 
the Is(2)/Is(o) ratio is given by the para excitation tem- 
perature Te^: 
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where A^ is the Einstein A coefficient for sponta- 
neous emission of level i (Ag(2) — 2.75 xlO~^, As(o) — 
2.94x10-", As(i) = 4.76xlO-Q |Wolniewicz et all [i998'). 
Ui is the frequency of transition, gi is the statistical 
weight of the level (gs(2)=9, gs(o)=5): and AE is the 
energy difference between the emitting levels (= 1172 
K). We plot the J=4 level excitation temperature so de- 
rived in Figure 6 (marked as S(2)/S(0)). The excita- 
tion temperature is remarkably uniform over the disk at 
Te:^(2,0)=204±8.5K. 

Between the para and ortho species, the line intensity 
ratios have the additional factor of the relative abun- 
dance, the ortho to para ratio, o/p. For a high temper- 
ature (T > 200 K) gas in local thermodynamic equilib- 
rium, this ratio is given by the ratio of the statistical 
weights, of the two species, o/p = 3. Assuming the o/p 
ratio is 3, then averaged over the eleven positions for 
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which the S(2) hne is detected (hereafter, the "S(2) re- 
gion"), the average extinction corrected S(1)/S(0) and 
S(2)/S(l) hne intensity ratios (1.92 and 0.59) yield level 
excitation temperatures of Tea;(l,0) = 125 K, Tea:(2,0) = 
204 K, and Tea;(2,l) = 388 K respectively. The level ex- 
citation temperatures are monotonically increasing with 
J, which is permitted if there is more than one emit- 
ting components within the beam, each with a differ- 
ent gas excitation temperature. For temperatures less 
than ^ 200 K , the cquihbrium o/p ratio is less than 3. 
iBurton et al.l fl992) calculate the equilibrium o/p ratio 
as a function of gas temperature. Making use of this 
work, we arrive at a self consistent solution for Tex (1,0) 
= 134 K with o/p = 2.30. Since the S(1)/S(0) line inten- 
sity ratios are fairly uniform across the galaxy, we assume 
this o/p ratio applies everywhere, and use it to compute 
Te3;(l,0). Notice that T^xl'i-fi) is not affected by the o/p 
ratio, and Te2;(2,l) is much greater than 200 K, so that 
o/p ^ 3. Figure 6 shows the computed excitation tem- 
peratures across the disk of the galaxy. As we found for 
Tea;(2,0), Tea;(l,0) is also remarkably uniform across the 
plane with a mean value of Te2:(l,0)=132±6K. Tea;(2,l) 
shows more variation, but this is largely not statistically 
significant. Excluding the data point 6 kpc to the NE 
of the nucleus, the mean value is Tea:(2,l) = 372 ±41 K. 
The position at 6 kpc NE of the nucleus is statistically 
higher excitation, with Te2,(2,l) = 556 -I-127/-87 K. 

4.2.2. Warm H2 Gas Mass Fraction: Enhancement of the 
Northern Ring 

Since the line emission is optically thin, the line in- 
tensity yields the gas column of the emitting level, N„: 
I = hi^s(n) Ag(n) Nu/(47r), whcrc As(n) is the Einstein 
A coefficient for the S(n) transition at frequency ;^s(n) 
for n = 0,1,2. Averaged over the eleven positions from 
which the S(2) line is detected, the column densities are 
N2 = 1.22 X 1020cm-2, N3 = 9.02 x lO^^ cm-^, and 
N4 = 6.81 X 10"^^ cm~^. The total column density of 
molecular hydrogen is given by: N^oi = (N„/g„) Z(T) 
g-(£;„/(fcT)) ^j^gj.g g^ ^ (o/p)/(2J+l), and Z(T) is the 
partition function: 
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The symbols have their meanings as before. Since the 
level excitations are not all the same, there must be more 
than one warm molecular gas component. The simplest 
model has two gas components, a "warm" component, 
which will emit most of the observed S(0) line, and a 
"hot" component that will emit most of the observed 
S(2) line. Both components will emit in the S(l) line. 
With a two component model there can be no unique so- 
lution as there are 4 variables (2 components, each with 
a Tea; and N(H2)), and just 3 observables (the lines). 
We model the line emission averaged within the "S(2) 
region" . The minimum mass solution will have a high 
gas temperature for the warm component. The largest 
gas temperature permitted by the data is the observed 
1-0 excitation temperature. For this temperature, all 
of the observed S(0) and S(l) line emission arises from 
the warm component. Therefore, there can be no self- 
consistent solution for all three H2 lines since the warm 



component is too cool to emit much in the S(2) line, and 
hot component, which will emit most of the the S(2) line, 
will also emit in the S(l) line. Using the partition func- 
tion, and the relationship between line intensity and col- 
umn density (Ij,j_i = hz^j_j_iAj^j_iN,7/(47r)), we have 
explored the gas temperature, column density parameter 
space for self consistent solutions, and find the minimum 
mass solution has gas temperatures of T^arm — 127 K 
(o/p ratio shifts to 2.2), and Thot = 1700 K, and H2 
column densities of N(H2, ujarm) = 3.8 xlO^^ cm~^, and 
N(H2,/iot) = 8.2 xlQis cm"2. Within this model, 75% 
and 25% of the observed S(l) line emission arises from 
the warm and hot components respectively. The warm 
component emits less than 5% of the observed S(2) line, 
while the hot component emits less than 1% of the S(0) 
line. The case where 25% and 75% of the observed S(l) 
line emission arises from the warm and hot components 
respectively is obtained for gas temperatures of T^arm — 
107 K (o/p ratio shifts to 1.65), and Thot = 465 K, and 
H2 column densities of ^(ii2,warm) = 5.8 xlO^^ cm"'^, 
and N(H2. hot) = 3.4 xlO^® cm~^. For this case, less than 
1% of the observed S(2) (S(0) line) arise from the warm 
(hot) component. For all cases, the warm component 
completely dominates (more than 99%) the H2 column 
density. 

We use the minimum mass self consistent model (75/25 
split of S(l) line flux into warm and hot components) to 
calculate the total column density in H2 line emitting 
molecular gas (esssentially the warm component from 
above) across the disk of NGC891. The results, along 
with the column traced in CO(l-O) are plotted in Figure 
7a. The warm molecular gas column peaks at the nucleus 
at a value of ~ 5.26 x 10^^ cm~^, and is above 2 x 10^^ 
cm~^ out to ± 8 kpc from the nucleus. Figure 7b shows 
the ratio of the warm molecular gas to that traced in its 
CO(l-O) line emission (within the Spitzer beam), i.e. the 
fraction of the total molecular gas that is emitting in the 
H2 lines. This fraction is modest ('^13%) and nearly con- 
stant in the inner galaxy from 3 kpc NE of the nucleus 
to 6 kpc SE of the nucleus, but rises to ~ 30% in the 
outer galaxy. This could be a true enhancement of the 
warm molecular gas mass in the outer galaxy, or it could 
be due to a change in the conversion factor, X between 
CO intensity and H2 column density with distance from 
the nucleus. If the X factor were larger by a factor of 
two in the outer galaxy compared with the inner galaxy, 
at first glance, this could explain our apparent change in 
the warm/cold molecular gas fraction. However, we do 
not believe this is the case, as there is a significant NE- 
SW asymmetry in the warm/cold molecular gas fraction 
apparent within the molecular ring (r = 3 to 7 kpc). Av- 
eraged over the ring positions, the NE portion of the ring 
has twice the warm molecular gas fraction (~ 28%) as 
the SW portion of the ring (~14%). The total molecular 
gas as traced by CO is very clo se to the same in both 
sides of the ring (jScoville et al.r i993): it is the warm 
molecular gas column that is enhanced (Figure 7, top). 
To cancel out this effect, one would need to invoke an 
X factor twice as large for the NE ring as for the SW 
ring. This enhanced warm molecular mass fraction is 
more likely the result of enhanced starformation in the 
NE ring, which is reflected in the stronger [CII] and [01] 
line emission observed to the NE as well (see Figure 8). 



Integrated over the galaxy, the warm molecular gas totals 
^ 6.1 X 10® M0, or about 16% of the total molecular gas 
mass as traced by CO (3.9 x 10^ M0) within our Spitzer 
beams (which extend ±5.5" of the galactic plane). 

4.3. Gas Heating 

The power to heat large amounts of molecular gas to 
temperatures above a hundred degrees can be provided 
by far-UV (stellar) photons, by the dissipation of molec- 
ular gas turbulent kinetic energy in cloud-cloud collisions 
and on smaller scales within clouds, by X-ray photons, or 
by cosmic rays. For normal star-forming galaxies, only 
the first two may provide enough power to be important. 

4.3.1. Photodissociation Regions - Motivation 

Figure 8 provides insights into the source for the heat- 
ing of the molecular gas. In this figure we plot the 
"cold" molecular gas mass (traced in its CO emission) 
the "warm" molecular gas mass (traced in its II2 line 
emission), and the [CII] 158 /zm and the [01] 63 fine- 
stru cture lines obtained with the ISO LWS spectrome- 
ter (|Brauher. Dale fc and Helou Il2008( ). We also include 
weak detections of the [01] 146 /im at 6 positions along 
the galactic plane that we obtained through close inspec- 
tion of the LWS L02 scans of the line obtained within 
the ISO d ata archive. These data were c alibra ted as de- 
scribed in IBrauher. Dale fc and HelouH (|2008[ ) so as to 
be consistent with their [CII] and [01] 63 /im fluxes, and 
their [Nil] 122 /i fluxes which we use below. The cold and 
warm molecular gas distributions have been smoothed to 
the significantly coarser spatial resolution (^75") of the 
ISO LWS beam for proper comparisons, and the trac- 
ers are normalized to their nuclear values (excepting the 
[01] lines which are normalized to their peaks) for ease 
of comparision. The cold molecular gas mass is centered 
on the nucleus with a significantly narrow distribution 
along the galactic plane than the warm molecular gas 
mass. In addition, the warm molecular gas mass distri- 
bution is asymetrical, with near nuclear values extending 
to the molecular ring, 4.6 kpc to the NE of the nucleus. 

The [CII] line predominantly arises from photodissoci- 
ation regions formed on the surfaces of molecular clouds 
exposed to the far-UV (6eV < hv < 13.6 eV) radiation 
from nearby OB stars or the general interstellar radia- 
tion field. The [CII] line is also an important coolant for 
low de nsity ionized gas, and "atomic clouds", but most 
(> 70 % [Stacev et al.lll99lllAbel et al.ll2005t lOberst et al.l 
l2006f ) of the [CII] flux from external galaxies likely arises 
from PDRs, so tha t the line is an ind icator of OB star 
formation activity (jStacev et al.lll99ll ). The very close 
correspondance between the [CII] distribution and that 
of the warm molecular gas is strongly suggestive of a 
PDR origin for the II2 line emission. 

The ionization potential for O (13.62 eV) is essentially 
the same as that of H, so that the FIR [01] 63 and 146 /im 
line emission also arises from neutral gas within PDRs. 
The [01] lines have significantly higher critical densities 
(~ 4.7 X 10^ and 9.4 x 10"* cm'^ for the 63 and 146 /im 
lines respectively) than that of the [CII] line (~ 2.8 x 
10^ cm~^) so that the [01] lines trace denser PDRs. The 
strong peaking of the [01] 63 /zm line line on the north- 
ern molecular ring is strongly suggestive of enhanced gas 
densities in this region. 



The remarkable spatial correlation between the warm 
II2 mass and the [CII] line including the extension to the 
NE suggests a common origin in PDRs. Indeed, warm 
molecular gas is expected within PDRs. The penetra- 
tion of carbon ionizing photons into neutral gas clouds 
is limited by the extinction of these photons by dust to 
visual extinctions. Ay '^ 3 mag (corresponding to an 
Nh ~ 6 X 10^^ cm~^), and along h alf this column, hydro- 
gen i s typically molecular in form (jTielens fc HoUenbachl 
Il985l) . 

4.3.2. Photodissociation Regions: Constraints on Go 

We have argued for a PDR orgin for much of 
the H2 line emission. As such, one might expect 
a correlation between the observed H2 line emission 
and that of an independent tracer of PDRs, the 
PAH 7.7 ^m emission feature (e.g. iHabart et al.l 
(|200llPeeters. Spoon, fc Tielens I ((2001 . In Figure 9, 
we compare the PAH 7.7/im emission with the H2 line 
emission, by plotting the H2-to-PAH emission ratio along 
the plane of NGC891. The PA H emission is meas ured 
from IRAC 8 /im imaging data (jWhalev et al.ll2009D us- 
ing the same apertures than those used for the H2 spec- 
troscopy. To compute the PAH emission from the 8 
/im flux, we use the same effective bandwith {\v = 
13.9 X 10^2 Hz ) as iRoussel eFall (|200l . The H2 line 
emission is the sum of the S(0), S(l) and S(2) lines. All 
tracers were corrected for extinction (see section 3.3). 
The extinction at 7.7 /im is similar to that at 12.3 and 17 
/im (|Drainell2003D . so that the H2 to PAH emission ratio 
is not very sensitive to the sp e cific c hoice of the extinc- 
tion c urve. iRigopoulou et al.l ([2002) and IRoussel et al.l 
(|2007f ) have shown that the H2 emission of star forming 
galaxies is tightly correlated with the PAH emission. Fig- 
ure 9 extends this correlation to the distribution of both 
tracers within a galaxy. The H2 to PAH emission ratio 
is remarkably constant within the NGC 891 disk, with a 
mean value only slightly larger than that measured for 
the SINGS galaxies. 

A striking outcome of the SINGS data analysis is that 
the H2/PAH emission ratio is independent of the dust 
heating as measured by the 24 /im to FIR dust emis- 
sion ratio. IRoussel et al.l (|2007D argue that this result is 
easiest to understand if the H2 and dust emission come 
both from PDRs, and are powered by a common energy 
source, the UV light from young stars. However, from 
their discussion it is not fully clear how this interpre- 
tation fits wit h the modeling of th e dust spectral energy 
distribution in lDraine et al.l ( 20071 ). where the PAH emis- 
sion comes mainly from interstellar matter in a low UV 
radiation field (typically G < 10). Here, we quantify the 
constraint set b y the H2 /PAH emis sion ratio on the UV 
field. Using the IDraine fc Li I (|2007D models available on 
line, we computed that the PAH emission, as it appears 
in Figure 9, is 
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where the dus-to-gas mass ratio is assumed to be 
Galacti c and qpAH is the fraction of the dust mass in 
PAHs praine fcLill200l . The default v alue of qpAH 
corresponds to the median value derived by (jDraine et al.l 
,2007. 1 for metallic galaxies in the SINGS sample. The 



model output may be combined with the observed value 
of H2/PAH {^ 0.012), to estimate the effective radiation 
field, G, within the warm H2 gas. For NGC 891, we 
find L(H2)/M(H2,warm) - 1/20, so that putting this ah 
together we have: 
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Here we have left in an explicit scaling factor for the 
warm molecular gas mass. Clearly, this calculation only 
provides a rough upper limit for G, since it assumes that 
the PAH emission arises exclusively from PDRs contain- 
ing the warm H2 gas, while some PAH emission is ex- 
pected from the H I and cold H2 gas. However, within 
the current understanding of PAH emission, the H2 to 
PAH comparison therefore suggests that the bulk of the 
warm H2 mass (the S(0) emitting gas) is located in PDRs 
exposed to modest far-UV radiation fields. We use this 
as a constraint in our PDR models that we develop be- 
low. It is important to note that this constraint on G 
only applies for the gas that radiates in the S(0) and 
S(l) lines, not the gas radiating in the S(2) line. The 
S(2) line emission could arise from a much higher G en- 
vironment without increasing the PAH flux very much 
since the mass of this component will be so small. 

4.3.3. Photodissociation Regions: Modeling 

iKaufman et al.l (|2006[ ) have made available a grid of 
PDR model outputs, which include the expected rota- 
tional line emission from H2 as a function of G, the 
strength of the local far-UV {& <hv < 13.6 eV) radiation 
field normalized to the Habing field (Gg = 1.6 x 10~^ erg 
s~^ cm~^), and the cloud number density n. The model 
computes the local heating and cooling, chemistry and 
radiative transfer as a function of depth from the cloud 
surface into its dust-shielded core. Diagnostic line ratios 
and line intensities appropriate for a single face-on slab 
geometry are posted on-line, and available interactively 
in the "Photodissociation Region Toolbox (PDRTJ3'- 
Table 2 lists our observed H2 S(0), S(l) and S(2) lumi- 
nosities and those of the [CII] , and [01] lines integrated 
over the galaxy. We estimate that the errors for the H2 
luminosities are roughly given by their corrections for 
extinction, and high latitude flux missed by the Spitzer 
beam. Figure 10 (left) plots our observed galactic aver- 
age H2 line ratios as a function of G and n from PDRT. 
The model results do not give a good fit to the NGC 891 
line ratios. The best fit is for high density (n> 6 x 10^ 
cm~^) and moderate far-UV radiation fields (G ^ 200 to 
400). 

Figure 10 (center) plots the [01], [CII] and ([01] 63 /.tm 
-I- [CII]) /FIR continuum ratios as averaged over the 
galaxy onto a PDR diagnostic diagram from PDRT. For 
this FIR data set we make corrections to the observed 
luminosities of our tracers. The first correction is for 
optical depth and geometry. The PDRT models are ap- 
propriate for face-on single slab geometry. Since the ISO 
beams encompass regions ~ 2.5 kpc in size at the dis- 
tance to NGC 891, there is an ensemble of clouds within 
our beam. For a spherical cloud, the projected surface 

^^ See dedicated web site at |http://dustem. astro. umd.edu] 



area is only 1/4 the total surface area, so that a beam 
filled with spherical clouds will have 4 times the inten- 
sity as a unit filling factor single slab if the lines are 
optically thin. If all the tracers are optically thin, taking 
the line ratios eliminates this geometric factor. Although 
it is almost certainly true that all the H2 lines are op- 
tically thin, optical depth may be an issue for the FIR 
lines. In fact, the 63 imh [01] is most li kely optically 
thick (jStacev et al.lll983l: iTielens fc Hollen bach 1985) in 
star formation regions, so that the observed [01] line in- 
tensity should be corrected for geometry, multiplying by 
as much as a factor of four for high optical depth in the 
line in spherical geometry. For our study, we make a 
more modest correction, multiplying the 63 fiui line by 
a factor of two. We take the [01] error boundary to be 
50% due to this factor of two correction, and note that 
this error bound will encompass geometry effects up to 
a factor of 3. The [CII] and [01] 146 /.tm lines are likely 
optically thin, or thinish as is the FIR continuum, so we 
make no optical depth correction to their fluxes. 

Next, we correct the observed [CII] line emission for 
that fraction expected to arise from diffuse ionized gas. 
This correction is perhaps best done by comparing the 
observed [CII] line flux to that of the ground state fine 
structure lines of N+ at 205 and 122 /xm. With an ion- 
ization potential of 14.5 eV, these lines only arise from 
within ionized gas regions, and one can show that the 
[CII] /[Nil] line ratios are a measure of the observed frac- 
tion of the [CII] li ne that arises from the io nized inter- 
stellar medium (cf. IPetuchowski fc Bennettl[T993 : Heilei 
ri99l lObCTst et al.l |2006D. We use the [Nil] 122 ^m lines 
obser ved by the ISO LWS (jBrauher. Dale fc and HeloiTI 
|2008() , averaged over the plane of the galaxy. The [Nil] 
122/xm line to [CII] line ratio in NGC 891 is ~1:9: some- 
what smaller than that of the Milky Way Galaxy (1:6, 
iWright et al.lll99ll ). Assuming Galactic C/N abundance 
ratios the line ratio indicates that between 50% (low den- 
sity limit) and 4% (high density li mit) of the observe d 
[CII] comes from ionized gas regions (l Oberst et al.ll2006l) . 
Here we will assume 27% of the [CII] line arises from HII 
regions as found for t he low density (ue'^ 30 cm~^) gas 
in the Carina nebula ()Oberst et al.l 120061) . We take the 
error bound for the [CII] luminosity to be equal to the 
fraction (27%) we have subtracted off. 

Finally, not all of the FIR continuum luminosity from 
galaxies will arise from dense PDRs. A significant frac - 
tion can arise from the d iffuse ISM. lAlton et al.l ()1998f ) 
and iPopescu et all (|2004[) have fit "warm" (PDR) and 
"cool" (diffuse) fractions of the total FIR luminosity from 
NGC 891. The warm dust likely arises from H II regions 
and PDRs, while the cool dust arises from the cold intcri- 
ors of molecular clouds and the diffuse ISM (Draine et al] 
I2007D . We therefore take --50% of the total FIR luminos- 
ity of NGC 891 to arise from PDRs. Finally, we correct 
L_F/fl to Lao-ioooAim which is appropriate for the PDRT 
models (see footnote to Table 2). Considering these cor- 
rections, we take the error bound for the FIR luminos- 
ity to be 50%. Putting the FIR line and continum cor- 
rections together, the PDR cooling line/cooling contin- 
uum ratio integrated over NGC 891 is ([01] + [CII])/FIR 
~1.7% (Table 2). 

Let us now examine the overall properties of NGC 891, 
integrated along the galactic plane. Taken within them- 



selves, the S(2)/S(0) and S(1)/S(0) line ratios (1.08 and 
1.69 respectively) indicate G -200-400, and n >6 x 10'* 
cm-3. However, the [01] 146Aini/[CII], [01] 63/im/[CII] 
and ([OI]-t-[CII])/FIR ratios indicate lower far-UV fields 
with best fit near G— 100, and a significantly lower den- 
sity with best fit near n— 3x10^ (Figure 9, center). The 
solution spaces for the H2 and FIR line models do not 
overlap. The closest approach for the two models is ob- 
tained for similar Gs (-- 100 to 200), but the FIR hues 
require a lower gas density (u'^ 10"* to 4 xlO'^ cm~^, as 
G ranges from 100 to 200) than the II2 lines (n> 6 x lO'' 
cm~^). If we imagine the PDRs to be the surfaces of 
externally illuminated molecular clouds, then since the 
H2 lines arise from deeper within the PDR than the FIR 
lines, it is likely they will arise from denser gas regions, 
so perhaps this is an acceptable solution. If so, then the 
beam filling factors of the two components should match 
as well. For the FIR component, the average [CII] line 
intensity (less 27% for ionized gas) referred to an 11" re- 
gion along the plane of the galaxy is 2.8 x 10"'* erg s~^ 
cm~^ sr~^. At the closest approach solution between 
the H2 and FIR models (n- 4 x 10^ cm'^ and G-200), 
the [CII] intensity predicted by PDRT is - 2.2 x lO"* 
erg s~^ cm~^ sr~* so that the beam filling factor for the 
FIR lines is 0fir linos --l-S. For n- 6 x 10"*, G~200, 
PDRT predicts a S(0) line intensity -- 1.0 x 10"^ erg 
s~^ cm^^ sr^^. The extinction corrected value within 
our 11" beam averaged over the plane of the galaxy is 
— 1.7 X 10~^ erg s~* cm~^ sr~^, so that the beam filling 
factor is ^Ua uncs ~ 1-7, in fair agreement with the FIR 
value. 

4.3.4. Is some of the H2 ermsston powered by the dissipation 

of turbulence ? 

The H2 lines, the FIR lines, and the FIR continuum 
are consistent with PDR models but only if there is a 
rather large (factor of —10) density gradient between the 
(denser) H2 emitting region and the (less dense) FIR 
line emitting regions. In addition, the requisite far-UV 
field strength is significantly larger than the value we 
derive using the observed PAH distribution in section 
4.3.2 above. These difficulties motivate us to investi- 
gate the possibility that a fraction of the observed H2 
line emission arises from the dissipation of turbulence 
in molecular clouds. This possibility has been proposed 
to ac count for diffuse H2 line emission in the Milky 
Way (' Falgarone et al.l j2005") and the galaxy-wide shock 
in Stephan's Quintet (JGuiUard et al. 2009 ). 

The S(1)/S(0) line ratio is fully consistent with the far- 
UV fields and densities as derived from the FIR lines and 
continuum (see Figure 10, left). What drives the H2 so- 
lution is the S(2) line, which is a facto r of six too bright. 
The C-s hock models of Dra i ne et al.l (Il983fl and, more 
recentlv iFlower fc Pineau-Des-ForetsI (|2010l ) are in good 
agreement, and indicate that the S(2) line is brighter 
than the S(l) line for Vshnck > 15 km s~^. F or example, 
using the iFlower fc Pineau-Des-Foreti (|2010l ) model, for 
shock velocities, Vghock ~20 to 30 km s~^, with preshock 
density n^ = n(H)+2n(H2)= 2x10'* cm'^, the S(2) line 
is ~ 1.6 x as bright as the S(l) line, and 75 times brighter 
than the S(0) line. If 80% of the observed S(2) line flux 
were to arise from such C-shocks, and the remainder from 
PDRs, then 30% of the S(l) line would also arise from 



shocks, and the PDR line ratios would shift to S(2)/S(l) 
~0.17, and S(2)/S(0)- 0.22. The allowed values for G 
and n for this PDR plus shocks solution are shaded blue 
in Figure 10 (right). Clearly there is now a large region 
of overlap between the PDR solution for the FIR lines, 
and that for the H2 lines. This solution space is centered 
at about G~60 ± 40, n~ 5 +5/-3 xlO^ cm'^ (Figure 
10, center). The revised filling factors are reasonably 
self-consistent with t/fpiR lines ~ 1-9, and 0h2 linos ^ 2.8. 
With Vshock ~ 25 km s~*, there is significant cooling in 
other molecular hydrogen lines: the S(3) and S(4) lines 
are predicted to be 4~7 and 4 times brighter respectively 
than the S(2) line. The total luminosity in the molecu- 
lar hydrogen lines from the shocked gas would be very 
large: Ln^ ^ 10^ Lq. For a 25 km s~^ C-shock, the 
common cooling lines of [01] and CO are more than 100 
and 15 times less important energetically than those of 
H2 , and H2O is predicted to have a total cooling ^ 
that of H2 ([Flower fc Pineau-Des-Foretsll2010t ). A simi- 
lar analyis of the region centered on the north-east ring 
has a solution centered on modestly enhanced excitation, 
Gne ring ^^80, Ti^E ring ~ 8 X 10'^ cm"'^, and rcquircs 
somewhat higher area filling factors (j)-piK iinos ~ 2.7, and 
0H2 linos ~ 3.6 to account for the enhanced [01] 63 /im 
line emission there. 

Leaving aside the question of the H2 excitation, we can 
assess the contribution of mechanical energy dissipation 
to the H2 emission, from the point of view of energet- 
ics. Is the dissipation rate of the interstellar medium 
turbulent kinetic large enough to provide a significant 
contribution to the H2 emission? Various groups have 
discussed the dissipation of turbulence using numerical 
magneto-hydrq dynamic sim ulations ( Stone et al. 1993; 
iMac Low! 1 19991 e.g.). As in iBradford et a l. (2003), one 
can rearrange equation (7) from IMac Low (1999) to pa- 
rameterize the energy released by the dissipation of tur- 
bulence per unit mass in terms of the dispersion of the 
turbulent velocity, v^ms , and the scale length of energy 
injection, A^: 



L 

— = 2.5 X 

M 



10-3 f 



VlOkms- 



3 /lOOpc 



— (5) 

Mo ^ ^ 



The default values for Yrms and A^; correspond to the 
upper end of the line widt h-size relation for gian t molec- 
ular clouds in the Galaxy (jSolomon et al.|[l987l ). We in- 
voke shocks to account for 80% of the S(2) luminosity, 
and 25% of the S(l) line luminosity. By dividing the sum 
of these line luminosiites 9.5 x 10^ L©, by the molecular 
mass 4.4 x 10^ M©, derived from the CO(l-O) luminos- 
ity, we obtain LH2/MH2 = 2.2 x 10"^^. Therefore, the 
expected turbulent dissipation rate is greater than the lu- 
minosites of the lines, so that a significant contribution 
from turbulent energy dissipation is likely. 

In summary, the observed H2 line emission from NGC 
891 is well fit by a model that combines a low excitation 
PDR (G ^ 60, n~5x 10'^ cm^^) ^ffiWy dissipation of tur- 
bulent energy within molecular clouds. The (adjusted) 
[01] (63 and 146 ^m) and [CII] 158 ^m fine structure 
line fluxes, (half) of the observed FIR continuum fiux, 
the H2 S(0) line fiux and 70% of the S(l) line flux arise 
from the PDR component, while the H2 S(2) line and 
30% of the S(l) line flux arise from gas heated by turbu- 
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lent dissipation. It is interesting to note that based on 
observations of the H2 l ines along a hne of sig ht within 
the Milky Way galaxy, iFalgarone et alj (|2005[ ) reached 
much the same conclusion: the S(0) line emission may 
be accounted for by UV heated H2 in the diffuse inter- 
stellar medium and low excitation PDRs, but the UV 
emission cannot be the unique energy source, because 
their PDR models fall short of the observed S(l) and 
S(2) lines by large factors (~2-10). As we outline here, 
they propose that the S(l) and S(2) H2 line emission 
is powered by the localized dissipation of turbulent ki- 
netic energy, albeit in a low density gas (uh ^ 50cm~^) 
within the diffuse interstellar medium. The presence of 
a diffuse component of H2 line emission has been con- 
firmed by additional Galactic observations ([Hewitt et al.l 
[20091 : [Goldsmith et al. 1 1201(11 ) . 

It should be clear that the solutions presented here are 
by no means unique. On galactic scales there will cer- 
tainly be a range of PDRs along the line of sight with 
varying degrees of excitation, that could mimic the ob- 
served line ratios, perhaps negating the need to invoke 
shocks. Furthermore, a single parameter shock model is 
also simplistic: there will very likely be a superposition 
of shocks of various parameters only any given line of 
sight. What we present here is intended to be only a 
representative solution. 

5. CONCLUSIONS 

We have mapped the distribution of the S(0), S(l), 
and S(2) pure rotational lines of H2 along the plane of 
the edge-on galaxy NGC891. The lines are remarkably 
bright. After making an extinction correction, the to- 
tal line flux from the three transitions is 2.82 x 10^ L© 
or ^0.09% of the TIR luminosity of the galaxy. Unlike 
the ISO results that suggest a change in the S(1)/S(0) 
line intensity ratio in the outer regions of the galaxy, we 
find the line ratios are nearly constant along the plane of 
the galaxy indicating uniform excitation of the gas. In 
particular, there is no reason to invoke a massive cold 
component of the molecular gas in the outer galaxy as 
the ISO results appear to require. Through analytical 
modeling, we find the bulk of the mass traced in its II2 
line emission has an excitation temperature near 125 K. 
This "warm" molecular component is a substantial frac- 
tion of the total molecular gas mass: about 10 to 15% in 
the inner (3 kpc NE < r < 6 kpc SW) galaxy rising to 
more than 30% in the outer regions. Integrated over the 

galaxy, Mwarm H2 ^ 16% Mcold CO- 

We compare our II2 observations with detailed PDR 
models to show that the observed line intensities and 



ratios are consistent with high density (n> 6 x 10"* 
cm~^) molecular clouds exposed to moderate far-UV 
fields (G-200-400). We refine our PDR modeling by 
using available ISO [01] and [CII] data, which indicate 
smaller far-UV fields (0^^5-300), and lower gas densities 
(n^ 2 X 10^ — 10^ cm~^). The solutions spaces for the 
II2 and FIR line PDRs do not overlap, but a reasonably 
close solution is obtained for G~ 200, and a large (factor 
of ten) density gradient between the FIR emitting gas 
and the H2 emitting gas. This solution may be possi- 
ble if the clouds are externally illuminated, but the high 
field invoked by this solution is in contradiction to the 
more modest far-UV fields (G~ 50) one would expect by 
comparing the observed H2 line luminosity to that of the 
7.7 /im PAH features. A consistent solution is obtained 
within a two component model. The first component 
is a low excitation PDR, with G ^ 60, and the second 
component is heated by the dissipation of turbulent ki- 
netic energy. We find a very good fit for the H2 and FIR 
lines and continuum if all the observed S(0) line emis- 
sion arises from the PDR component, and 80% of the 
S(2) line arises from the turbulent component. The S(l) 
line emission is split 70% into the PDR component, and 
30% into the turbulent component. This solution is only 
representative of a range of solutions that could be ob- 
tained by allowing more PDR and/or shock components 
with different excitation parameters. 

This is the first analysis of the H2 rotational line emis- 
sion from a galaxy together with the FIR lines and con- 
tinuum in the context of PDR models. A complication 
with this modeling is the relatively large beams and mod- 
est sensitivies of the ISO spectrometers for detection of 
the FIR lines. These issues will become much less im- 
portant in the near future with the arrival of FIR data 
from the 3 m class Stratospheric Observatory for Infrared 
Astronomy (SOFIA) and the Herschel Space Telescope 
spectrometers and photometers. 
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TABLE 1 
H2 Line measurements along NGC 891 



ID 


RA 


Doc 


Distance^ 


S(2) A12.3Atm 


S(l) A17.0Atm 


S(0) A28.2^m 




J2000 


J2000 


arcmin 


xlQ-" (Wm-2) 


xlQ-" (Wm-2) 


xlO-17 (Wm-2) 


26 


2:22:44.10 


+42:26:2.450 


5.65 


<1.40 


<1.24 


<1.44 


25 


2:22:43.35 


+42:25:33.85 


5.53 


<1.47 


<1.39 


<2.10 


24 


2:22:42.37 


+42:25:06.67 


4.97 


<1.74 


<1.37 


0.77+0.26 


23 


2:22:41.44 


+42:24:40.09 


4.42 


<1.67 


<1.58 


2.41+0.63 


22 


2:22:40.78 


+42:24:07.65 


3.87 


0.67+0.11 


1.45+0.11 


5.08+0.62 


21 


2:22:39.83 


+42:23:39.90 


3.32 


<1.89 


1.77+0.30 


5.77+0.85 


20 


2:22:38.38 


+42:23:13.93 


2.76 


1.09+0.09 


2.49+0.16 


8.84+1.07 


19 


2:22:37.59 


+42:22:46.22 


2.21 


4.11+0.45 


4.10+0.48 


10.91+0.68 


18 


2:22:36.47 


+42:22:18.77 


1.66 


2.15+0.20 


4.13+0.68 


11.37+0.45 


17 


2:22:35.47 


+42:21:50.29 


1.11 


2.79+0.63 


3.82+0.50 


11.07+0.48 


16 


2:22:34.67 


+42:21:21.43 


0.55 


2.14+0.59 


3.88+0.57 


10.44+0.80 


1 


2:22:33.41 


+42:20:56.88 


0.00 


3.43+0.19 


6.86+0.57 


14.64+0.58 


2 


2:22:32.37 


+42:20:25.85 


-0.55 


<1.58 


4.35+0.59 


11.66+0.54 


3 


2:22:31.02 


+42:19:59.96 


-1.11 


3.23+0.61 


5.15+0.30 


13.65+0.51 


4 


2:22:30.16 


+42:19:31.63 


-1.66 


1.90+0.50 


3.80+0.40 


9.29+0.59 


5 


2:22:29.15 


+42:19:03.34 


-2.21 


1.82+0.55 


2.51+0.75 


7.39+0.44 


6 


2:22:28.07 


+42:18:35.81 


-2.76 


1.28+0.11 


2.81+0.21 


6.82+0.68 


7 


2:22:26.99 


+42:18:08.74 


-3.32 


<2.70 


2.61+0.44 


7.20+0.34 


8 


2:22:25.98 


+42:17:41.93 


-3.87 


<2.19 


0.92+0.24 


3.45+0.64 


9 


2:22:24.88 


+42:17:14.69 


-4.42 


<2.40 


1.58+0.50 


4.24+0.53 


10 


2:22:23.74 


+42:16:46.54 


-4.97 


<1.41 


0.56+0.08 


1.86+0.21 


11 


2:22:22.70 


+42:16:19.70 


-5.53 


<1.64 


<1.05 


<2.24 


12 


2:22:21.65 


+42:15:52.55 


-6.08 


<1.44 


<1.07 


<1.66 



Offset of each pointing along the disk of NGC 891, from the pointing on the nucleus (position ID — 1) 



TABLE 2 
Integrated Line Luminosities in NGC 891 



Line 



Luminosity 

(L0) 



Corrected Luminosity* 
(L©) 



Error 



Ha S(0) 
Ha S(l) 
Ha S(2) 
[OI] 63 fim 
[OI] 146 Mm 
[CII] 158 iim 



5.44x1015 
8.94x10'' 
5.26x10*5 
4.67x10'^ 
4.74x10'^ 
1.40x10* 
1.51x10^" 



7.18x106 
1.32x10'' 



7.80x10' 
9.34x10'' 
4.74x10*5 
1.02x10' 
1.13x10^0 



24% 
24% 
32% 
50% 
50% 
37% 
30% 



^ H2 correction factor from extinction (1.16, 1.32, and 1.32 for the S(0), S(l) and S(2) lines respectively), and high latitude emission missed by our 
beam (1.16 correction factor). The [OI] 63 fini corrected by factor of 2 for geometry and [CII] line is corrected by 0.7 for HII region contribution. 
The FIR (42.5 to 122.5 fim) luminosity is divided by two to account for cirrus emission, then multiplied by 1.5 for flux longward of 122.5 /J.m as 
required by the PDRT models. The factor of 1.5 is appropriate for galaxies in the Revised Bright Galaxy Survey l lSanders et al.|[2603i) 
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Fig. 1. — A Spitzer/IRAC 8/im image of the NGC891. The IRS pointings are indicated with circles of 25" 
where at least one H2 line was detected and listed in Table 1 are marked with their corresponding number. 
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Fig. 2. — Zoomed regions near the H2 lines of the final IRS spectra at each one at the pointings listed in Table 1, where at least one of 
the S(0), S{1), and S(3) lines were detected. The spectra are arranged so the left, center, and right columns refer to the S(2), S(l) and 
S(0) lines respectively, and each row refers to a position on the sky (as listed in Table 1) 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Offset from nucleus (kpc) 

Fig. 3. — a) Observed Spitzer H2 line intensities as a function of distance from the nucleus. Positive offsets are to the north-east, b) 
Observed Spitzer/ISO S(0), and S(l) line flux ratios normalized to the ISO beam sizes. 
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Fig. 4. — a) The column densities of hydrogen nuclei contained in H2 (derived from CO(l-O), see section 3.3), HI, and the sum of the 
two as a function of offset (in kpc) along the galactic plane from the nucleus, b): Extinction correction, C(A), for the observed molecular 
hydrogen line emission. Note that the correction is the same for the H2S(1) and H2S(2) lines. 
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Offset from nucleus (kpc) 

Fig. 5. — a) Extinction corrected Spitzcr H2 line intensities as a function of distance from the nucleus, b) S(1)/S(0), S{2)/S(0) and 
S(2)/S(l) line intensity ratios calculated from the extinction corrected data in panel a. 
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Fig. 6. — Level excitation temperatures determined from the extinction corrected line intensity ratios. For Texi^fi) we assume the o/p 
ratio is given by thermodynamic equilibrium at T = 135 K: o/p = 2.3. 
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Fig. 7. — a) Column densities in H2 as traced by the CO(l-O) line (dotted line, "cold" H2) and as traced by the S(0) and S(l) lines (solid 
line, "warm" H2). Error bounds are in grey, b) The ratio of "warm" to "cold" H2 column densities. We have enclosed the regions of the 
NE and SE ring by dotted lines. 
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Fig. 8.— Cuts in the [CII] (1 58 fim), and [OI] (63 and 146 fim) lines along the plane of NGC 891 (from the ISO satellite, 
IBrauher. Dale fc and Helou I I I2008I )), superposed on our derived cuts in the "warm" and "cold" H2 gas which have been smoothed to 
the spatial resolution of the ISO beam (75"). There is excellent correspondence between the warm H2 and the [CII] line including the 
asymmetry to the north of the nucleus strongly suggesting a PDR origin for the H2 line emission. 
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Fig. 9. — Ratio between the H2 line emission, summed over the S(0), S(l) and S(2) Unes and the PAH 7.7/im feature emission (measured 
from the IRAC 8/^m imaging data as in Roussel et al. (2007)) along the plane of NGC891. The ratio is remarkably constant though the 
inner parts of the galaxy and rises by about a factor 2 at the profile edges. The range of values measured in SINGS star forming galaxies 
(mean value and la scatter) is shown for comparison. 



Fig. 10. — PDR solutions for the G and n from the Photodissciation Region Toolshead (PDRT) | |http://dustem. astro. umd.ed u/ll. (left) 
Solutions for the H2 lines. The red shaded region indicates allowed solutions when all of the line emission is trom fUtis. (center) Solutions 
for the FIR lines and the ([OI] + [CII])/FIR continuum luminosity ratios. The blue shaded area indicates the allowed solutions with the 
various corrects to the tracers described in the text. The red shaded area is the solution for the H2 lines from the PDR/shock model on 
the right, (right) Solutions for the H2 lines after subtracting off a shock contribution. The red shaded region indicates allowed solutions 
when 80% of the S(2) line emission, and 30% of the S(l) line emission is assumed to arise from shocks. 



